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SUMMARY 

Round trip missions to Mars have been investigated 
to define representative launch windows and associated 
AV requirements. The 1982 inbound and the 1986 out- 
bound Venus swingby missions were selected for analy- 
sis and serve to demonstrate the influence of the char- 
acteristics of the heliocentric trajectories on the launch 
window velocity requirements. This report presents 
results indicating the effects on the launch windows of 
velocity capability, transfer technique, and of the incli- 
nation, eccentricity, and insertion direction of the orbit. 
The analysis assumed a circular parking orbit at Earth 
and considers both circular and elliptical parking orbits 
at Mars. Use of one-, two-, and three-impulse trans- 
fers were investigated. The three-impulse transfer 
employs an intermediate elliptic orbit of 0.9 eccen- 
tricity, For all cases, insertion at planet arrival was 
into an orbit coplanar with the arrival asymptote and 
any required plane change was performed during the 
planet departure phase. 

The minimum AV requirement to transfer from a 
circular parking orbit to a hyperbolic asymptote occurs 
when the orbits are coplanar and the maneuver is per- 
formed at periapsis of the hyperbola. The study indi- 
cates that, using a three-impulse transfer, the AV 
penalty for non-coplanar departures, is no more than 
5-10% above the minimum coplanar requirements. 
Therefore, use in mission analysis of the coplanar AV 
requirements would not result in large errors if three- 
impulse transfers are acceptable. Use of fewer im- 
pulses significantly increases the error. Similar 
characteristics occur for elliptical parking orbits. How- 
ever, due to the low coplanar AV's, they provide a 
longer launch window for a given total AV capability. 


In this paper, a comparison of three different tech- 
niques for performing the turning maneuver from cir- 
cular orbit during launch windows at both Earth and 
Mars is made. The techniques consist of optimum one- 
and two- impulse transfers, and a restricted three- 
impulse transfer utilizing an intermediate elliptic orbit 
from which the required plane change is made. In 
addition, an assessment of using three similar transfers 
from elliptical parking orbits at Mars is made and the 
results are compared to the circular orbit results. 

To be most meaningful, the departure analysis 
should be mission oriented. Thus, the orbit escape 
techniques were applied to two specific Mars stopover 
missions. These two missions encompass the four 
trajectory classes of primary interest: i.e., direct in- 
bound and outbound, and Venus swingby inbound and out- 
bound transfers. Therefore, while the quantitative re- 
sults are not directly applicable to other Mars missions, 
the conclusions can be generalized to other round trip 
trajectory modes by utilizing combinations of these 
transfer trajectories. 

The paper is divided into three major sections. The 
first discusses the orbit transfer techniques which were 
utilized and provides parametric data in support of the 
assumptions and approximations made in the analysis. 
The second section provides the criteria for the mission 
analysis and defines the characteristics of the 1982 in- 
bound Venus swingby and of the 1986 outbound Venus 
swingby which were chosen as representative Mars 
missions. The coupling of the missions with the trans- 
fer techniques is contained in the third section, which 
provides contour maps of constant velocity increments 
as functions of inclination and staytime in orbit. 


INTRODUCTION 

Preliminary analysis of interplanetary missions is 
generally performed with the assumption that the day- 
to-day launch window penalties can be closely approxi- 
mated by the change in the coplanar departure velocity 
requirements. Since a parking orbit at Earth and/or 
the objective planet will generally be used, a unique 
orbit inclination will exist during the launch period of 
interest. The time dependence of the departure hyper- 
bolic asymptote and the regression of the orbit plane 
due to planetary oblateness allow essentially coplanar 
departures for only short time periods. In general, a 
non-zero angle will exist between the asymptote direc- 
tion and the orbit plane, and the velocity increment 
necessary to accomplish the required turning maneuver 
can be of such magnitude as to render the coplanar 
assumption invalid. — 
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METHOD OF ANALYSIS 

Orbit Geometry 

In considering the launch window problem for a 
given mission, it is first necessary to establish the 
relative positions of the orbit plane and the hyperbolic 
escape asymptote with time. For the analysis of plane- 
tary launch windows, it was assumed that the parking 
orbit at planet arrival was coplanar with the arrival 
asymptote. The resulting orbit elements are shown in 
Figure 1. A planet- centered right-hand coordinate 
system with Z axis at the north pole and X axis at the 
planet vernal equinox was chosen as indicated. The 
arrival asymptote or hyperbolic excess velocity vector 
is defined conventionally as a planet centered vector 
with right ascension, p, and declination, 6. The orbit 




elements of interest, that is, the inclination, i, and 
longitude of the ascending node, Q, are related by 


sin (p - Q) 


tan 6 
tan i 


magnitude of (p - Q) is greater than 90°. This departure 
is called northern injection. The other plane has the 
spacecraft departing below or south of the excess velocity 
vector. This case has a (p - ft) less than 90° and is a 
southern injection. 


for all i > 6, while the argument of periapsis, co, is 
measured from the ascending node. 


It can be seen that for each inclination there are two 
orbit planes coplanar with the arrival asymptote. These 
two orbits can be distinguished by the relative positions 
of the spacecraft approach vector at the time of the 
final midcourse maneuver and the planet centered 
excess velocity vector. If the spacecraft approaches 
above the excess velocity vector, it moves initially 
toward the north pole of the planet and the magnitude 
of (p - ft) will be greater than 90°. This orientation 
shall be referred to as a northern insertion. In the 
other case, the spacecraft approaches below the excess 
velocity vector (i.e., toward the south pole) and the 
magnitude of (p - SI) is less than 90°. Figure 1 illus- 
trates this configuration which is called a southern 
insertion. 


The situation at some time after arrival is illus- 
trated by Figure 2. The perturbation due to planet 
oblateness causes the orbit plane to regress about the 
planet in the manner shown. That is, the inclination of 
the orbit remains unchanged, the longitude of the ascend- 
ing node changes by A £2, and the argument of periapsis 
changes by Acu where, to first order in the planet 
oblateness, 1 
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P = a(l - e 2 ) 

At that time , the orbit plane makes an angle 1^ with the 
departure asymptote. 2 

sin loo = cos i [sin 5^ - tan 6 r cos 6^ cos (p r - p d - fit)] 
±[sin 2 i - sin 2 sec 6 r cos 6 ^ sin (p r - p d - fit) 


It is this angle which must be compensated for during 
the departure maneuver. 


The orbit geometry at Earth was analyzed in a 
similar manner except that the initial orbit (corre- 
sponding to the orbit at arrival for the planetary case) 
was chosen in such a manner as to be, through regres- 
sion, coplanar with the departure asymptote at the nomi- 
nal departure time; that is, zero plane change was 
chosen for nominal departure. As with the planetary 
orbits, two orbital planes exist which satisfy this con- 
straint. These planes are distinguished by the relative 
position of the spacecraft departure and the departure 
excess velocity vector at the nominal departure date. 

If the departure is above the excess vector, then the 
spacecraft departure is toward the north pole and the 


Impulsive Analysis 

The assumption was made for this analysis that the 
velocity increments are applied impulsively. The first 
effect of this assumption is the neglect of gravity losses. 
Gravity losses, however, are relatively small for these 
missions and can easily be approximated, or even ne- 
glected without significant error. 

It can be shown that the orientation between the 
parking orbit and the escape hyperbola is preserved 
whether the velocity is added impulsively or through a 
finite thrusting time. This is indicated by the example 
in Figure 3 where the velocity increment required to 
escape to a hyperbolic excess speed of 6.0 km/sec from 
Mars is shown as a function of the turning angle, that is, 
the angle between the line of apsides and the departure 
asymptote. This example is for a departure from an 
orbit with an eccentricity of 0.9 and 1000 km periapsis 
altitude. The initial thrust-to- weight ratio for the 
finite thrusting maneuver was optimized as a function 
of the true anomaly at the start of thrusting, and varied 
from about 0.03 near apoapsis, to 0.4 near periapsis. 

The primary difference between the results for 
impulsive thrusting and for finite thrusting is in the 
position on the parking orbit at which the velocity is 
added. For finite thrust, the start of thrusting occurs 
prior to (i.e., leads) the true anomaly for impulsive 
velocity addition. That is, to obtain the same excess 
speed and direction as in the impulsive case, the finite 
thrusting must be approximately centered about the 
position of the impulsive velocity addition. Typical lead 
angles are shown in Figure 4. Here, the velocity incre- 
ment to escape to an excess speed of 6.0 km/sec from 
Mars is shown as a function of true anomaly around an 
orbit with an eccentricity of 0.9 and a periapsis altitude 
of 1000 km. The true anomaly shown for the finite thrust 
is that at the start of thrusting. It can be seen that for 
equal velocity increments (and hence equal turning 
angles) the start of finite thrusting must lead the posi- 
tion for impulsive velocity addition by approximately 
5° to 15° in true anomaly for all true anomalies except 
near apoapsis. Near apoapsis, the associated slow 
rotation and finite thrust result in a near zero lead 
angle requirement. 

In summary, the assumption of impulsive velocity 
increments is attractive due to the simplification it 
permits and appears justified due to the small resultant 
errors. 

Circular Parking Orbit 

A comparison of three methods of transfer from a 
circular parking orbit to a hyperbolic asymptote is 
provided. These methods are optimal one- and two- 
impulse direct transfers and a restricted three- impulse 
transfer using an intermediate elliptical orbit. The 
techniques are briefly delineated in this section. 

One-Impulse Transfers - The single- impulse 
transfer problem is one of solving for the minimum 
velocity increment (AV) to achieve the desired 


hyperbolic asymptote from a specific orbit. The exact 
solution for minimum AV results in an expression which 
does not permit a closed form solution. 2 ^ However, in 
Reference 3 an approximate solution is developed which 
provides very good agreement with the exact solution 
for interplanetary velocity ranges. 

Two-Impulse Transfers - The two-impulse trans- 
fer technique considered in this study utilizes one im- 
pulse at departure from the circular parking orbit and 
the second impulse at an "infinite" distance from the 
planet, i.e., at the sphere- of- influence. The solution 
used was developed in Reference 3 and employs an opti- 
mal distribution of velocity change and angle change 
between the two impulses. Reference 2 considered a 
transfer with an optimum distribution of velocity be- 
tween the two impulses and with the second impulse sup- 
plying the entire plane change. The parametric results 
are nearly the same as the optimal transfer employed 
here. 

Figure 5 presents the parametric data for one- and 
two- impulse transfers. The figure indicates the regions 
of velocity and plane change requirements where each is 
superior in terms of minimum total AV requirement. 

Three-Impulse Technique - The three- impulse mode 
considered in this study utilizes an in-plane tangential 
impulse to insert the spacecraft from a circular parking 
orbit into an intermediate highly- elliptical orbit followed 
by a plane change at apoapsis to rotate the orbit plane so 
as to be coplanar with the departure asymptote. The 
third impulse is then provided tangentially at the appro- 
priate position on the intermediate ellipse so as to place 
the spacecraft on the required escape hyperbola. The 
position of the first impulse on the circular parking orbit 
is varied parametrically and from these results the 
optimum position is determined by a numerical search. 
In general, the optimum position is such that the plane 
change maneuver is minimized without an inordinate 
penalty for an off-periapsis departure from the inter- 
mediate ellipse. A complete analysis of this three- 
impulse technique has been performed. 4 

A brief study was made of off-apoapsis plane-change 
maneuvers on the intermediate orbit with the motive of 
rotating the line of apsides in order to reduce off- 
periapsis departure penalties. It was found that for 
highly eccentric orbits very little was to be gained by 
such maneuvers and that the optimum plane-change 
position was very near apoapsis. This result is pri- 
marily a reflection of the low off-periapsis departure 
penalties for highly elliptic orbits. The off-apoapsis 
plane- change option was not retained further in this 
study. 

The effect of the eccentricity of the intermediate 
ellipse upon the total departure velocity requirements 
is shown in Figures 6(a) and 6(b) for Earth and Mars, 
respectively. These data are for a typical hyperbolic 
excess speed of 6 km/sec and for an orbit altitude of 
1000 km. Curves for various values of the angle be- 
tween the orbit plane and the departure asymptote, 1^, 
are shown. It can be seen, especially for large out-of- 
plane angles, that high eccentricities are very attractive. 

Parametric data for the three- impulse departure 
requirements as a function of hyperbolic excess speed 
are shown in Figures 7(a) and 7(b) for Earth and Mars, 


respectively. These data are for an initial circular 
orbit altitude of 1000 km and an intermediate orbit 
with an eccentricity of 0.9. Again, various values of 
the out-of-plane angle, 1^, are shown. The minimum 
in departure velocity requirements with hyperbolic 
excess speed which occurs for the Earth case at 
rather low values of excess speed as exhibited by the 
loo = 30° line is due to an increasing off-periapsis de- 
parture penalty as excess speed is decreased. That is, 
for a given 1^ (non-zero), there is a minimum excess 
speed below which a periapsis departure is not possible. 

Elliptical Parking Orbits 

For the elliptical parking orbit, the sum of the 
arrival and the departure AV requirements must be 
minimized. This approach is necessary since the posi- 
tion of periapsis significantly affects both of the AV 
requirements and the coplanar arrival maneuver will 
establish that position. In general, both the arrival and 
the departure will be performed off-periapsis to achieve 
the minimum total AV. For all transfer techniques, the 
orbit insertion point was varied from +90° to -90° true 
anomaly on the parking orbit. 

As with the circular orbits, three transfer tech- 
niques were considered. 5 They are an optimum single 
impulse and constrained two impulse and three impulse 
transfers. The techniques are briefly described below. 

One-Impulse Technique - In this technique, since 
both the departure hyperbolic excess vector and the 
characteristics of the parking orbit at departure are 
fixed, the non-coplanar, non-tangential departure 
maneuver true anomaly is iterated to establish the 
departure point of minimum AV requirement for each 
arrival true anomaly. In general the departure is 
off-periapsis and has a non- minimum plane change 
angle. 

Two-Impulse Technique - The first impulse is 
applied at apoapsis to perform the required plane 
change. The second impulse is applied tangentially 
near periapsis and provides the required velocity 
change to achieve the desired excess velocity. Thus, 
this technique is similar to the last two impulses of the 
three- impulse transfer from circular orbit. 

Three-Impulse Technique - The three-impulse 
transfer for elliptical orbits is similar to that for cir- 
cular orbits. The first impulse is applied at periapsis 
to increase the orbit eccentricity; again, the value of 
0.9 was used for the transfer ellipse. The second im- 
pulse performs the necessary plane change maneuver 
at apoapsis and the third impulse applied near periapsis 
provides the additional velocity increment necessary to 
achieve the desired hyperbolic trajectory. 


MISSION SELECTION CRITERIA 

Data to aid in selecting reasonable mission pro- 
files is contained in Ref. 6, which presents mission 
characteristics for both direct and Venus swingby 
stopover missions to Mars for each Earth-Mars oppo- 
sition period from 1980 to the year 2000. 

The results of that study indicate that with the 
exception of the 1984 and 1997 oppositions the swingby 






mode is more attractive than the direct mode. Propul- 
sive velocity requirements are consistently lower, in 
some cases affording reductions of thirty percent. 

Earth entry speeds are also consistently lower, afford- 
ing reductions of up to fifty percent. These benefits 
are realized for mission durations of about 500 days, 
which represents an increase of about 20 percent over 
the durations for direct missions (based on thirty-day 
stopovers). These mission durations were obtained by 
minimizing the product of propulsive velocity require- 
ments and mission duration. These "nominal 11 missions 
from Ref. 6 have mission durations about 25 percent 
shorter than minimum energy missions at the expense 
of less than ten percent increase in total propulsive 
velocity requirements. As will be demonstrated later, 
such a mission selection procedure contributes to low 
sensitivity of the mission parameters to launch delays. 
For the above reasons, and to encompass typical in- 
bound and outbound swingbys, the launch window analy- 
ses were conducted for the "nominal" missions from 
Ref. 6 for both the 1982 inbound swingby and the 1986 
outbound swingby. 

Since the intent of this paper is to assess the in- 
fluence of launch delays on mission characteristics, 
one approach to the selection of the round trip trajec- 
tories could have been the specification of a nominal 
stopover time. This would have necessitated a further 
specification on the parking orbit characteristics at 
Mars to ensure that some other criterion be met, e.g,, 
that the nominal departure from Mars be coplanar, or 
that the plane-change velocity penalty be minimized 
throughout some arbitrary staytime, etc. The approach 
taken here i3 to vary both stopover time and orbit in- 
clination parametrically and to divorce the selection of 
the outbound leg from the selection of the inbound leg. 

Representative Earth departure and Mars arrival 
conditions must nevertheless be specified. For both 
mission years considered, the outbound leg chosen is 
that of the nominal mission employing a 30-day stop- 
over. The Mars departure date was then varied to 
reflect launch delays. While this may seem restrictive, 
it is in fact quite reasonable. For outbound swingbys, 
regardless of mission duration or stopover time, varia- 
tions of no more than about 10 days arise in either the 
optimum Earth departure or Mars arrival dates. For 
inbound swingbys, the variation in Mars arrival date 
can be on the order of 50 days. However, while the use 
of the optimum arrival date would reduce the Mars 
coplanar departure velocity requirements, it does not 
significantly influence the plane- change requirements 
for the various stopover times. 

This paper is not directly concerned with system 
characteristics, but recognizing that Mars arrival 
velocities should not be unconstrained, the outbound 
legs are selected so that the sums of the Earth depar- 
ture and Mars arrival velocity increments are equal to 
the nominal values insofar as possible. When this con- 
stant value can no longer be maintained for the longer 
staytimes, the trajectories are then selected on the 
basis of minimizing the velocity sum. Note that this 
procedure is possible by virtue of the manner by which 
the nominal missions were defined. If the nominal 
missions had been selected on the basis of minimum 
energy, then all delays could conceivably cause the 
velocity requirements to increase. Figures 8(a) and 
8(b) illustrate the variation of outbound leg parameters 
for the 1982 and 1986 missions, respectively. 


The technique of maintaining constant velocity is 
also applied successfully to the definition of the inbound 
legs in the face of launch delays at Mars as shown in 
Figures 9(a) and 9(b). In this case, only the Mars de- 
parture velocity is considered. Constraints could also 
be imposed on Earth entry speeds, of course, but as can 
be seen they remain quite low. To maintain either a 
constant Mars departure velocity or minimum Mars 
departure velocity for staytimes in excess of about 
60 days, it is necessary for the inbound leg heliocentric 
transfer angle to exceed 180°. Thus, the inbound leg 
duration must increase accordingly. 


NUMERICAL ANALYSIS 

In order to provide an assessment of the perfor- 
mance of one-, two-, and three- impulse transfers, the 
techniques described for transfer analysis were applied 
to the two selected missions. This section presents, 

(1) the results of analysis of the Earth launch window 
which employs only circular parking orbits, and (2) the 
results of the analysis of Mars 1 launch windows em- 
ploying both circular and elliptical parking orbits. 

The format selected for presentation of the launch 
window data is a contour map, showing lines of constant 
propulsive velocity requirement on an orbit inclination 
versus departure date plot. Thus, for a given orbit, the 
launch window for a given propulsive capability is 
readily visible. These data are obtained by cross 
plotting the complete set of transfer data which is com- 
puted for fixed inclinations as a function of departure 
date. The data necessary to prepare the contour maps 
are included in References 4 and 5. 

Earth Launch Window 

The launch site considered for these missions was 
the Kennedy Space Center (KSC). Therefore, the orbit 
inclinations which can be obtained without a plane 
change are required to be between the site's latitude 
(« 28°) and the range safety constraint (about 50°). The 
Earth launch window analysis was, therefore, limited to 
posigrade orbits of 30°, 40°, and 50° inclination. 

1982 Inbound Swingby - The results for the 1982 
launch are shown in Figures 10(a) and 10(b) for a cir- 
cular orbit of 300 km altitude and using an intermediate 
eccentricity of 0.9 for the three-impulse transfers. 

Since the variation of the Earth departure hyper- 
bolic excess velocity for this opportunity is from 4 to 
5 km/sec, Figure 6 indicates that the two- impulse trans 
fer will have lower AV requirements than the one- 
impulse for plane-change angles (I*,) over 15°. Figure 10 
shows that the two impulse does indeed lower the AV 
requirement at a given departure date. However, for 
reasonable velocity penalties associated with the plane 
change (i.e., 10-20 percent of the minimum coplanar 
departure AV) the increase in launch window size 
(« 1 day) is not of sufficient magnitude to make the 
two- impulse transfer of practical interest. 

The use of three- impulse transfers is seen to sig- 
nificantly increase the available launch window. For a 
AV of 4.5 km/sec, the three- impulse window is over 
24 days for the southern injection (Fig. 10(a)). For the 
same AV, only a 2-day window is available using one 
impulse transfers. To achieve a 24-day window with the 
one-impulse transfer would require a AV of 6 to 7 km/sec . 


Thus, the increase in launch window flexibility using a 
three-impulse transfer is such as to make it worth con- 
sideration in mission analyses even with the complexity 
of making three separate maneuvers. 

The southern injection orbit rotates relative to the 
departure vectors such that two coplanar departures 
occur during the departure period studied. Departures 
between these points require plane changes and thus 
AV penalties. The three-impulse transfer significantly 
reduces this plane change AV penalty thus opening the 
region to low AV departures. Comparison of Fig- 
ure 10(a) for the southern injection and Figure 10(b) 
for the northern injection reveals the effect that the 
orbit direction can have upon the launch window. 

For the northern injection, Figure 10(b), the orbit 
rotation relative to the departure vector variation is 
such as to result in only one coplanar departure in a 
departure period studied. A second coplanar departure 
occurs at an earlier time. From Figure 8, this earlier 
departure has a higher Earth departure AV requirement. 
The plane change penalty for this orbit orientation is 
also higher than for the southern injection. Therefore, 
the use of a three- impulse transfer for the northern 
injection orbit results in AV requirements about 10% 
higher than for the southern injection orbit in 1982 for 
the same window length. 

1986 Outbound Swingby - The Earth departure data 
for the 1986 opportunity is shown in Figure 11(a) for 
the southern injection orbit, and 11(b) for the northern 
injection. Conclusions drawn above as to the relative 
value of one, two, and three-impulse transfers for the 
1982 opportunity apply here as well. The effect of orbit 
orientation is different, however. For this opportunity, 
it is the northern injection orbit, Figure 11(b), whose 
rotation, relative to the departure vectors, allows two 
coplanar departures over the time period studied. 

Again, use of the three-impulse transfer opens the 
entire departure region between these coplanar points 
to low AV requirements. 

Mars Circular Orbit Launch Window 

The minimum inclination which will allow coplanar 
arrival and departure is defined by the magnitude of the 
larger of the declinations at arrival and departure. This 
inclination is 28.7 and 20.8 degrees for the 1982 and 1986 
missions, respectively. The maximum inclination is 
180° minus the minimum inclination, i.e., a retrograde 
orbit. The inclinations were, therefore, varied from 
30 to 150 degrees for both missions in 10 degree inter- 
vals, assuming a circular parking orbit of 1000 km 
altitude. 

Only data for one- and three- impulse transfers are 
shown in the following sections, since the reduction in 
AV achieved through use of the two- impulse transfer, 
when a reduction occurred, was small compared to the 
total AV requirement. The three-impulse transfer 
technique again employs an intermediate orbit with an 
eccentricity of 0.9 to perform the plane-change 
maneuver. 

Use of the circular parking orbits allows the arrival 
maneuver to be performed tangentially at periapsis of 
the arrival hyperbola. The insertion AV is, therefore, 
the same for all staytimes, inclinations and insertion 
directions. This AV is 2.7 and 4.3 km/sec for the 1982 
and 1986 Mars arrivals, respectively. 


1982 Inbound Swingby - Figures 12(a) through (d) 
present the contour maps of the AV requirements for 
the 1982 mission. The figures show the data for both 
southern insertions and northern insertions as previ- 
ously defined. 

Figure 12(a) shows that for a one-impulse escape 
from a southern insertion orbit, three regions exist 
having a nearly coplanar departure AV requirement 
(5 km/sec). These are near-polar orbits for staytimes 
up to 50 days and both low inclination posigrade and 
high inclination retrograde orbits at 50 days staytime, 
plus or minus 10 days. In order to provide continuous 
departure capability up to about 80 days for all avail- 
able inclinations, a AV of 9 km/sec is required. The 
low AV region near 90 degrees inclination occurs since 
the arrival and departure right ascensions are nearly 
equal for about 60 days staytime. 

Three-impulse transfers from a southern insertion 
orbit, Figure 12(b), make available low inclinations for 
short staytimes (5 days) and a low AV (5 km/sec) capa- 
bility. When the AV capability is increased to 6 km/sec 
all inclinations obtainable can be achieved for staytimes 
up to 70 days. This is a significant improvement over 
the one- impulse transfer. A similar increase in avail- 
able staytimes does not occur by increasing the AV 
another 1 km/sec since the coplanar velocity require- 
ment rises rapidly after 60 days staytime (Figure 9(a)). 

For this 1982 mission, use of a northern insertion 
at arrival increases the launch window for a low to 
moderate AV capability. This is shown in Figures 12(c) 
and (d). The different insertion direction effectively 
provides an additional 20 days of low AV requirement 
as can be seen by comparing Figures 12(a) and (c). 

This orbital configuration allows a AV of 6 km/sec 
with a one-impulse transfer to provide staytimes for 
all inclinations of up to 40 days and for some inclina- 
tions of up to 70 days. The three- impulse transfer 
provides basically the same launch window capability 
for either orbit insertion since the plane- change penalty 
is inherently small. 

1986 Outbound Swingby - The data for the 1986 Mars 
launch windows are presented in Figures 13(a) through 
(d). Between 60 and 70 days staytime a discontinuity 
occurs in the return leg characteristics. This results 
from the trajectory selection since the central angle of 
the transfer conic crosses the 180° ridge. Solutions 
can, of course, be found in this region by utilizing a 
finer grid of stay times than was done here, or by con- 
straining the central angle. 

Figure 13(a) shows the data for a southern inser- 
tion with a single- impulse departure. The data does 
not display the symmetry shown by the 1982 mission 
since the arrival and departure right ascensions are 
at least 45° apart for the first 60 days of staytime. 

After 60 days, the symmetry is stronger. For near 
coplanar departure AV capability (3.0 km/sec) only 
small strips of the inclination- staytime map are avail- 
able. A AV of 4 km/sec opens up the low inclination 
region for staytimes up to about 80 days. However, to 
achieve polar orbits for this mission, AV’s of about 
5 km/sec would be required if only a one impulse 
transfer were used. 

Comparison of these southern insertion data to 
Figure 13(b) for a three-impulse transfer reveals the 
effect of reducing the plane-change penalty. With a 


three-impulse transfer, a AV of about 3.1 km/sec will 
allow use of any obtainable inclination for the staytime 
range shown. 

The effect of using a northern insertion at arrival 
for the 1986 mission is shown in Figures 13(c) and (d). 
For this mission, the additional low AV region at short 
staytimes, mentioned for 1982, is only about 10 days 
and does not significantly change the areas available 
for a given AV level. A AV of 5 km /sec is required 
for either insertion to provide a reasonably large 
inclination- staytime spectrum for single- impulse 
transfers. 

The effect of the different insertion upon the three- 
impulse data (Figure 12(d)) is negligible with a AV of 
about 3.1 km/sec required to open the entire region of 
inclinations and staytimes studied. 

Mars Elliptical Orbit Launch Window 

The analysis was performed for parking orbit 
eccentricities of 0.3, 0.5 and 0.7 for the 1982 and 1986 
miss ions. 5 Data for all three transfer techniques was 
developed, but only one- and three- impulse transfers 
for eccentricities of 0.3 and 0.7 are shown. In the case 
of the three- impulse transfer, the intermediate orbit 
for the plane-change maneuver has a 0.9 eccentricity. 

As with the circular parking orbits, the declinations of 
the arrival and departure vectors define inclination 
limits between which the possibility of coplanar arrival 
and departure exist. 

In the analysis of elliptical parking orbits two 
parameters, in addition to those for the circular parking 
orbit, are of importance. These are the orbit eccen- 
tricity and the location of the line of apsides relative to 
the arrival and departure asymptote. For circular or- 
bits, the arrival AV is a function of the interplanetary 
trajectory leg and orbit altitude only. However, for 
elliptical orbits, an off- peric enter insertion maneuver 
at arrival can reduce the off-pericenter requirement at 
departure and thus lower the sum of the AV’s at arrival 
and departure. The value of the optimal off-pericenter 
angle at arrival is a function of the inclination of the 
parking orbit and the desired staytime at the planet. 

In order to assess the launch window requirements, a 
nominal staytime is selected and the optimum insertion 
maneuver for that staytime is defined for each inclina- 
tion. That insertion maneuver is then held fixed for all 
staytimes to determine the available departure launch 
window. Re selection of the insertion maneuver would 
be expected to produce some gains, but based upon the 
lew cases which were investigated, changes which 
resulted in reduced AV's at long staytimes also re- 
sulted in increased AV's for short staytimes. In order 
to decrease the volume of data but still present the 
salient features of elliptical parking orbits, only data 
for the better insertion direction of each year as 
established by the circular orbit data is presented 
(i.e., northern insertion for 1982 and southern inser- 
tion for 1986), 

1982 Inbound Swingby - Figure 14(a) presents the 
total AV contours (arrival plus departure) for a one- 
impulse transfer for an elliptical orbit of 0.3 eccen- 
tricity as a function of inclination and staytime. If a 
coplanar periapsis arrival and departure were possible, 
a minimum AV of 6.6 km/sec could be achieved. As can 
be seen in Figure 14(a), two small regions exist for 


which the arrival and departure conditions are suffi- 
ciently near the absolute optimum that essentially no 
penalty occurs. A six percent increase (0.4 km/sec) 
in the total AV results in a staytime of at least 20 days 
for all inclinations. Using the total AV representative 
of minimum launch windows for circular orbits (7.7 km/ 
sec), a staytime of at least 40 days exists. This is to 
be compared to the 5 -day staytime for the one- impulse 
transfer from circular orbits. 

Data for a three-impulse transfer from a 0.3 eccen- 
tricity orbit are shown in Figure 14(b). As in the cir- 
cular orbit case, an increase in launch window length 
over the one-impulse transfer occurs. However, the 
increase is relatively less than for the circular orbits. 
This is due to the close positions of the arrival vector 
and the departure vector and their relative magnitudes. 
For this year, the arrival and departure pericenters 
are located so that both maneuvers can be performed 
with one impulse near pericenter since the plane-change 
requirements are small. A AV of 7.7 would produce a 
launch window of about 48 days for any inclination except 
near equatorial. 

Figure 15 presents the contour plot for a one- 
impulse transfer from a 0.7 eccentricity parking orbit. 
For this orbit, the minimum possible AV requirement 
is 5.45 km/sec which is nearly achieved in two regions 
of Figure 15. An increase of five percent in AV to 

5.7 km/sec allows a window of 50 days for near polar 
and high inclination posigrade orbits. The particular 
form of this AV contour in the retrograde region is of 
interest in that it indicates a departure window of up to 
50 days, but not for staytimes below 20 to 40 days at 
various inclinations. If this early depart restriction is 
of concern, then the retrograde AV requirements are 
higher than those for posigrade orbits. Use of a AV of 

7.7 km/sec results in a launch window of 75 days. 

The use of a three- impulse transfer from this orbit 
does not significantly change the launch window contour 
and is, therefore, not shown. 

1986 Outbound Swingby - Data for a one-impulse 
transfer from a 0.3 eccentricity parking orbit is shown 
in Figure 16(a). A minimum AV requirement of 6.0 km/ 
sec exists for the idealized coplanar, pericenter maneu- 
vers; a minimum of 6.1 was actually achieved. A AV of 
8.3 km/sec must be allowed before a generally available 
launch window exists for all staytimes. Even then a 
significant benefit over the one-impulse transfer from 
circular orbit does not exist. This occurs since the 
arrival and the departure asymptotes are over 45° dif- 
ferent in right ascension. As a result, the one- impulse 
case requires a serious compromising of the departure 
maneuver. 

Figure 16(b) presents the contour maps for the 
three-impulse transfer. In this case, as compared to 
1982, the plane change and the departure maneuvers 
having been optimized individually allows a significant 
increase in the available launch window. The use of a 
total AV representative of minimum launch windows for 
circular orbits (7.3 km/sec) results in a launch window 
of over 60 days for almost all inclinations. 

The contour map for a one-impulse transfer from a 
0.7 eccentricity orbit is shown in Figure 17(a), For this 
orbit, the minimum AV requirement would be 4.85 km/ 
sec. A value this low was not achieved anywhere on the 


f 


one-impulse transfer plot. A AV of 7.3 km/sec opens 
the region below 50° and above 120° for staytimes of 
100 days. A AV of 8.0 km/sec allows use of the entire 
contour map. 

The results of using a three- impulse transfer from 
the 0.7 eccentricity orbit are shown in Figure 17(b). 
Here a AV of 6.3 km/sec opens essentially the entire 
map to use, while a AV of only 5.3 km/sec would allow 
any retrograde orbit. 


CONCLUDING REMARKS 

This paper provides a summary of the variation in 
launch windows that can be achieved through use of 
circular and elliptical orbits and different transfer 
techniques. The data is presented for two opportunities 
which are representative of future round trip missions. 
The data can be used to trade-off AV requirements 
with the complexity of multiple-impulse transfers for 
all inclinations and staytimes of interest. 

To provide departure capability on any day for 
nominal staytimes up to 60 days at Mars and 20 days 
at Earth can result in AV penalties up to 50 percent of 
the nominal coplanar AV (arrival plus departure) if 
one- impulse transfers are employed from circular 
orbits. However, three-impulse transfers from 
circular orbits permits the same launch windows 
(i.e., 60 days at Mars and 20 days at Earth) for AV 
penalties on the order of 5-10 percent of the total AV 
requirement. 

Elliptical parking orbits do not lend themselves to 
such general conclusions due to the cross coupling of 
the arrival and departure maneuvers. Increasing 
eccentricity of elliptical orbits does, of course, reduce 
the total AV requirement for a given launch window. 

In addition, the use of three impulses to perform the 
transfer generally requires lower AV's than using a 
single impulse. However, orbit orientation can exist 
(e.g., 1982, e = 0.7) for which no gain is shown by the 
three-impulse maneuver. 

As pointed out under the mission selection criteria, 
supportable, but somewhat arbitrary, ground rules were 
established to define the interplanetary trajectories. 
Nominal orbit parameters were established based upon 
these criteria, thus determining the motion of the orbit 
as a function of staytime. With an orbit established, 
reselection of the interplanetary trajectory can be 
employed to vary the departure characteristics and 
possibly reduce the AV requirements for a particular 
window. Preliminary data in Reference 4 indicates 
that reselection of the direct legs does reduce AV 
requirements for single-impulse transfers. Little 
effect was apparent for three- impulse transfers due to 
their already low AV requirements. Reselection of the 
swingby leg was also ineffective since the swingby re- 
quirement severely restricted planet departure 
characteristics. 
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APPENDIX A 
NOTATION 

a semi-major axis, km 

AV propulsive velocity increment, km/sec 

e eccentricity 

H altitude, km 

1*3 minimum angle between departure vector and 

orbit plane, degs 

i inclination, degs 

J 2 second harmonic of planetary oblateness 

R planet radius, km 

t time, sec 

Voo hyperbolic excess velocity, km/sec 

x, y, z components of the planet centered coordinate 
system 

6 ^ declination, degs 

r unperturbed orbital period, sec 

p right ascension, degs 

9 longitude of ascending node, deg 

9 orbital regression rate, degs /day 

a; argument of periapsis, deg 

Subscripts 

c circular 

d departure data 

i intermediate ellipse for three-impulse 

transfer 

p periapsis 

r reference data 



FIGURE LEGENDS 


Figure 1 .- Orbit geometry at planet arrival. 

Figure 2 .- Orbit geometry at planet departure. 

Figure 3.- Effect of finite thrusting time on turning angle. 

Figure 4 .- Effect of finite thrusting time on ignition true anomaly. 

Figure 5 *- Optimum one-impulse and two- impulse transfers. 

Figure 6.- Effect of eccentricity on three-impulse transfers, (a) Earth 
orbit, (b) Mars orbit. 

Figure 7 *- Effect of excess speed on three-impulse transfers, (a) Earth 
orbit, (b) Mars orbit. 

Figure 8.- Outbound leg trajectory characteristics, (a) 1982 inbound 
Swingby. (b) 1986 outbound swingby. 

Figure 9 -~ Inbound leg trajectory characteristics, (a) 1982 inbound 
swingby (Mars arrival date = 244 5210 ). (b) 1986 outbound swingby 

(Mars arrival date = 244 65OO). 

Figure 10 .- Earth departure AV contours . (a) 1982 inbound swingby; 

southern injection, (b) 1982 inbound swingby; northern injection. 

Figure 11 .- Earth departure AV contours, (a) 1986 outbound swingby; 
southern injection, (b) 1986 outbound swingby; northern injection. 

Figure 12 .- Mars departure AV contours, (a) 1982 inbound swingby; 
southern insertion; circular parking orbit; one impulse, (b) 1982 
inbound swingby; southern insertion; circular parking orbit; three 
impulse, (c) 1982 inbound swingby; northern insertion; circular 
parking orbit; one impulse . (d) 1982 inbound swingby; northern 

insertion; circular parking orbit; three impulse. 

Figure 13 .- Mars departure AV contours, (a) 1986 outbound swingby; 
southern insertion; circular parking orbit; one impulse . (b) 1986 

outbound swingby; southern insertion; circular parking orbit; three 
impulse . (c) 1986 outbound swingby; northern insertion; circular 

parking orbit; one impulse, (d) 1986 outbound swingby; northern 
insertion; circular parking orbit; three impulse . 

Figure l 4 . - Mars AV contours, (a) 1982 inbound swingby; northern 
insertion; parking orbit eccentricity = 0.3; one impulse, (b) 1982 
inbound swingby; northern insertion; parking orbit eccentricity =0.3; 
three impulse. 

Figure 15 .- Mars AV contours; 1982 inbound insertion; northern insertion; 
parking orbit eccentricity = 0.7; one impulse. 


Figure l6.~ Mars AV contours, (a) I9S6 outbound swingby; southern 
insertion; parking orbit eccentricity = 0.3; one impulse, (b) 1986 
outbound swingby; southern insertion; parking orbit eccentricity = 0.3; 
three impulse. 

Figure 17 .- Mars AV contours, (a) 1986 outbound swingby; southern 
insertion; parking orbit eccentricity = 0.7; one impulse, (b) 1986 
inbound insertion; northern insertion; parking orbit eccentricity = 0.7; 
three impulse. 
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Figure l6.- Mars AV contours . (a) 1986 outbound swingby; southern 
insertion; parking orbit eccentricity = 0.3; one impulse, (b) 1986 
outbound swingby; southern insertion; parking orbit eccentricity =0.3; 
three impulse. 

Figure 17 .- Mars AV contours, (a) 1986 outbound swingby; southern 
insertion; parking orbit eccentricity = 0.7; one impulse, (b) 1986 
inbound insertion; northern insertion; parking orbit eccentricity = 0.7; 
three impulse. 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
AMES RESEARCH CENTER, MOFFETT FIELD, CALIFORNIA 



Equator, 


Departure 



Orbit plane 
at departure 

Orbit plane 
at arrival 


Figure 2. 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
AMES RESEARCH CENTER, MOFFETT FIELD, CALIFORNIA 



Mars departure 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
AMES RESEARCH CENTER, MOFFETT FIELD, CALIFORNIA 



K^/V, 


Figure 5. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
AMES RESEARCH CENTER, MOFFETT FIELD, CALIFORNIA 



16 



Vqo = 6 Km/sec 
H = IOOO km 


0 .2 .4 .6 .8 1.0 

Intermediate orbit eccentricity 


Figure 6(a) . 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
AMES RESEARCH CENTER, MOFFETT FIELD, CALIFORNIA 



Figure 6(b) . 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
AMES RESEARCH CENTER, MOFFETT FIELD, CALIFORNIA 


Total departure velocity, km/sec 



0 2 4 6 8 10 12 

Hyperbolic excess speed, km /sec 
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Figure 9(a) . 
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